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ABSTRACT 


The report describes the cne-~dimensional Lagrangian hydrodynamics computer 
progran, PUFF. The code is used primarily in the study of X-ray effects. 

In the past year it has been extensively revised and is now quite different 
from versions used outside the Air Force Weapons Laboratory. The major cal- 
culations in each subroutine are explained with a comrlete descrirtion of all 
input-output variables. Sample problems with the appropriate data deck are 
included to allow a user to become familiar with data arrangement ang to c ‘ck 
the program on his computer. 
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SECTION Tf 


INTRODUCTION 


The purpose of this report is to provide users of the PUFF code with an 
explanation of the latest version of the program. Severs] sample problems are 
included so that the user may become familiar with data input procedures and, 
at the same time, calibrate the program for his machine. The sample problems 
were designed to demonstrate a variety of input conditions and should be ex- 


amined only in that context, 


PUFF is a one-dimensional Lagrangian hydrodynamics program used primarily 
in the study of X-ray effects. The code uses finite difference methods to 
solve the equations of hydrodynamics using slab geometry and assuming one- 
dimensional flow through a l-centimeter-square pipe (reference 1). One may 
follow a pressure pulse which is propagated through 4a sample consisting cf up 
to six separate materials. A grid or mesh is superimposed on the materiais in 
the sample. The grid areas are called zones and the grid lines are called zone 
boundaries. The mesh is numbered, increasing from left to right. The pressure 
pulse normally starts from the left and progresses through the mesh to the 
right. When we speak of the front of the mesh, we are referring to the larger 
zone numbers. When we speak of the back of the mesh, we are referring to the 


smaller zone numbers. 


The index J is consistently used in the code to denote zone numbers. 
The materials in the sample are numbered consecutively from left to right. 


The index M is used when a material number is needed. 


The original PUFF was a derivative of a Livermore code (reference 2). 
Extensive work on PUFF was done at the Air Force Weapons Laborator; and ad- 
ditional work was done under contract by McAllister and Associates, Inc. 
During this development the code became quite large for its degree of sophis- 
tication. Many options as to types of problems run were incorporated and an 
extremely complicated and time-consuming edit routine was developed. Since 


this development took place over sa span of several years, many different 
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people made separate changes and no concentrated effort was made to delete 
options that were no loneer used. This made the program almost impossible to 


use for one not intimately familiar with it. 


In January of 1964 it was decided to make a version which would retain 
the essential elements of PUFF but eliminate all options which were not 
currently in use, and also to make the input-output as straightforward as 
possible. The monitor subroutine was eliminated and the edit routine was 
simplified by using the premise that it is more efficient to dump large amounts 
of required output on binary tape for later edit and plotting and to print only 


essential data during problem run time. 


Since the code was designed to be driven by one of two dissimilar methods, 
flyer plate (plate slap) and energy deposition, it was decided that a con- 
Siderable amount of computer time could be saved if the two options were 
handled in separate codes. The energy deposition code is now called PUFF and 
the plate slap code is called P PUFF. A great simplification in FORTRAN cod- 
ing resulted from this division as well as a reduction in core space and 
computer time required. During this revision cgs units were incorporated and 


the equation of state was completely divorced from the hydro loop. 


The present code consists of a main program where the hydrodynamics is 
accomplished and from which the other subroutines are called as needed. Eac’ 
routine has been extensively rewritten and simplified with the exception of 


rezone where only minor corrections have been made. 


PUFF and P PUFF are essentially the same except for the generator sub- 
routines. The section dealing with the generator subroutine points cut the 


major differences between the program. 


ae) 
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SECTION II 


MAIN PROGRAM 


Variables and counters used cnly in the main prugiam are initialized 
there instead of in the generator subroutine. Except for local variables 
teing initialized, all calculations in the main program are contained in the 
time or cycle loop. The time loop begins with the main hydro loop. This 
loop starts with the first zone, advances in time all quantities associated 
with that zone, and proceeds with each zone in sequence until all zones have 
been advanced to the new time. After the hydro locp is completed, numerous 
checks are made to determine if the problem is compiete, if a rezone is 
possible, or if one of the edit options is desired. if the problem is not 
complete, a new time step is calculated, the cycle counter N tis advanced 

1 


» and a naw cycle is started, 


The velocity U, the ordinate of the zone boundary X, the specific 
volume V, and the artificial viscesity Q are computed in a standard 
Lagrangian manner. The pressure P, energy E, and sound speed CS are cal- 
culated by use of an interpolation method using values returned from an equation- 
of~state subrountine. This method was developed at the Air Forte Weapons 
L.boratory and eliminates the requirement that the equation of state be in an 


analytic form. 


A tentative time step is computed using the Courant condition modified for 
the use of artificial viscosity. The actual time step used is the minimum of 
the tentative time step, 120 percent of the old time step, or 1 one-i.indredth 
of the energy deposit time. 


The Lagrangian edit records the time-pressure history of the first zone 
in the material to the immediate right (downstream) of the chosen zone boundary. 
This edit is accomplished in the main program since it is exercised every cycle. 
The zone designated as a JEDIT is not allowed to collapse in rezone. Ex- 
perience has shown it best not to designate consecutive zones as JEDITS unless 


one is very careful with the initial zone size. 
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SECTION IIlI 


GENERATOR 


bv Introduction 


The generator subroutine is called at the start of & problem to read all 
input data and initialize quantities needed in the problem. It is arranged 
in independent segments which are executed sequentially. All input values 


are printed out as a check of the data arrangement and for future reference. 


The zoning for both PUFF and P PUFF is computed using a sequence of geo- 
metric progressions. This is accomplished using two sets of variables RZ 
and NZ. The RZs are the common ratiss and the NZs are the zone indexes 


which terminate each progression. 


The only essential difference in the generators of PUFF and P PUFF is 
that PUFF determines an energy deposit rate SS while P PUFF sets a velocity 
profile. These items are discussed below. 


Qs Ener Deposition (PUFF 


The absorption coefficients are computed using the formula 


B 
uw = p'AA'(hey) - 


The values of AA and B are input numbers referenced to 1 kev that change 


with materiul and absorption edges. 


One may use the energy spectrum obtained fron a multiple black body source 
or may input en arbitrary spectrum. The arbitrary spectrum may be described 
by as many as 109 values of energy and corresponding h*v intervals. PUFF 
now uses 109 values of h‘v/KT in the black bedy calculations instead of the 
100 used in previous versions. This change was made to improve the high energy 
section of the curve. An energy deposit rate, 5S, is computed for each zone 


using the above values in the standard exponential absorption formula. 


2s Piate Slap (P PUFF) 
FP PUFF is designed to simulate the reactions caused by hitting a stationary 


plate with a piece of material (the flyer plate) moving at a given velocity 


ge eg Tey ne a 
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upon impact. The starting conditions for the problem simulate those at the 
instant of impact. 


The velocity of each zone in the flyer plate, excluding the last, is set 
equal to the input variable UZERO. Provision is made to smear the velocity 
drop, from UZERO to ZERO, over five zones. The velocity of the first zone in 
the target is set at an input fraction (UFACE) of the flyer velocity. UFACE 
is used to match the acoustical impedances of the flyer plate and the target. 


Zep 


Srp * 2p 


The velocity in the last zone of the flyer plate is the average of the flyer 


It may be computed from UFACE = > where 2 = * sound speed, 


Po 


plate velocity and the velocity of the first zone in the target. The velocity 
of the second zone in the target is one-half that of the first. A print of the 


calculated velocity profile is given. 
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SECTION IV 


REZONE 


ks Introduction 


At its best, rezone will introduce small perturbations in the solution. 
At its worst, rezone can corpletely destroy a problem. It is hoped that with 
judicious use the errors introduced will not affect the result by more than a 


few percent. The computer time saved is substantial. 


The rezone routine was designed for a single pulse moving through a mesh 
having initial zone size increasing with zonal index. Since available com- 
puter time limits a problem to a few hundred zones, rezone tries to distribute 
the zones where they are most needed. The routine attempts to keep the resolu- 
tion under the pulse constant by reducing zone size as the pressure peak 
approaches. Certain sections of the mesh are relatively inactive and will 
allow larger zone sizes. The collapse routine will combine two neighboring 
zones into one if certain conditions are met. The dividing and combining of 


zones is done in the divide and collapse routines respectively. 


[XN 
Wawsan 


lst zone ist zone 3rd zone zone of 20th 

in mesh left of to left max zone 

expanded max pres- of max pressure in front 

less than sure hav- pressure of zone 

20 ing less of max 
than 2/10th pressure 


max pressure 


Figure 1. Rezone Regions 
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2. Divide Routine 


The divide routine operates on twenty zones immediately in front of the 
zone of maximum pressure (region A, figure 1). To determine an appropriate 
Size for a test zone, the positive momentum of the main pulse and the maximum 
pressure are used to construct a theoretical triangular pulse which is assumed 
to be moving with velocity equal to the sound speed of the zone of maximum 
pressure, The pulse length of this theoretical triangular pulse is found and 
divided by an input number, ZPUL. This distance is used as the "optimum" 
zone size in the region. Any zone larger than this check size is divided in 
half. The divided zone is then checked to see if it should be redivided. 


3. Collapse Routine 


The collapse routine combines adjoining zones in regions B and C in 
figure 1. Region B is assumed to be more active than C. Region B extends 
from three zones to the left of the zone having maximum pressure to the first 
zone to the left having a pressure less than two-tenths maximum pressure. 
Region C begins at the left boundary of B and extends to the first zone 
in the mesh which is expanded less than 20 percent, or to the first zone if no 
zone is expanded 20 percent. In region B the optimum zone size computed in 


the divide routine is multiplied by one-half the relative volume (p)/o) of 


the zone of maximum pressure and this quantity is used for the optimum zone 

size in the regicn. Each zone in the region smaller than this check size is 
combined with an adjoining one. Region C uses a rather arbitrarily determined 
length for its zone size check. The length of the region is divided by the 
input variable JRZL and used as the optimum zone size. Each zone in the region 
smaller than this check size is combined with an adjoining zone. Zones desig- 
nated as material boundaries or JEDITS are not collapsed. The collapse routine 


is presently entered only once every 100 cycles. 
L. Comments 


If one is interested in the back part of the mesh (reflected shocks, etc.) 
at later times, JRZL should be increased to something like 100. A JRZL of 20 
will allow very large zcnes at late times which will effectively damp pulse 


movement. A ZPUL of 4&0 is adequate for most problems. 
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The divide routine may be bypassed in a problem by initializing ZPUL 
to zero. The collapse routine may be bypassed in a problem by initializing 


JRZL to zero. 
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SECTION V 


EDIT 


l. Introduction 


Elaborate edit routines usually consume a sizable percentage of the 
computer time used by a problem. This fact leads to a conflict between ob- 
taining sufficient edits to ensure adequate knowledge at all points of interest 
in the problem and the computer time that is available for the problem. Compli- 
cated edit routines also prove difficult in converting from one computer to 


another. 


The present edit routine is divided into two parts. The first part stores, 
at designated intervals during the problem, all zonal variables and other se- 
lected variables on binary tape for later processing. The second part of the 
routine prints a sufficient number of key variables to enable the user to 
determine during run time that the problem is progressing satisfactorily, and, 
after problem completion, to determine points where more detailed edits are 
wanted. This method has the advantage of using the minimum amount of computer 
time while providing the maximum amount of available data. It has the dis-~ 


advantage of requiring the user to wait longer for the final edits and plots. 


The Air Force Weapons Laboratory uses the Cal Comp automatic plotter for 
most PUFF data presentation. The most commonly used plot program simply plots 
the pressure versus distance for a given time, although more complicated 


programs exist. 


2. Output Variables 


All variables printed in the edit routine are defined below. 


N Cycle number. 

TIME Present problem time. 

DTN Present time step. 

JSTAR Zone number of the last zone in the mesh (from 


left to right) that has a velocity greater than 


1073 em/sec. 
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XJ PMAX 


JPMAX 


MV PULSE 


MV PREC 


MOMENTUM 


DTPP 


DTPULS 


ETOTAL 


EMVNEG 


EMVPOS 


L-BOUND 


The maximum pressure contained in the mesh. 


Coordinate of the zone boundary of the zone 


having maximum pressure. 


Zone number of the zone having the maximum 


pressure present in the mesh. 


Sum of the momentum of each zone from 
JPMAX+3 back to the first zone which as a 
negative velocity. 


Sum of the momentum of each zone from 
JPMAX+h out to JSTAR+1. 


Sum of MVPULSE and MVPREC. This expression 
is computed by averaging the zone mass of 
two adjcining zones and using the velocity 
of the boundary between. 


Pulse width computed using the expression 
EMVPP/PMAX. 


Pulse width computed using the expression 
EMVPL/PMAX. 


Total energy of the problem, computed by 
summing the kinetic and internal energy 


for each zone. The units are calories. 


Total negative momentum in the mesh. 
This expression is computed by averaging the 
velocities of two adjacent zone boundaries 


and using the zone mass of the zone between. 


Total positive momentum in the mesh. This 
expression is computed by averaging the 
velocities of two adjacent zcne boundaries 


and using the tone mass of the tone between. 


Coordinate of the left material boundary 
(X(1)). 


i 
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X(JBND1) 


X(JBND2) 


R-BOUND 


JFIN 


Coordinate of the boundary between the first 
and second material. It will be zero for one~ 


material problems. 


Coordinate of the boundary between the second 
and third material. It will be zero for two- 


material problems. 
Coordinate of the right material boundary. 


Zone number of the last zone in the problem. 


a1 
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SECTION VI 


EQUATION OF STATE 


a Introduction 


PUFF can now use any equation-of-state subroutine that will return 4 


pressure when presented an energy, density, and material index. 


The present equation of state is probably the greatest source of error 
in PUFF other than rezone. The pressure is computed using one of two functional 
equations: cone for solid regions and one for vapor regions. If the zone is 
compressed relative to ambient conditions, the solid state equation is used. 
If the zone is expanded relative to ambient conditions, the vapor state equa- 
tion is used. 


2. Solid Phase 


The solid phase equation is based on the Mie-Griineisen equation of state 
(reference 3). This equation appears in PUFF as 


3 


P= Cy + Due + Su + Gok , 


where u is p/egrl, G is the Criineisen parameter and E is energy density. 


C is determined by the product of ambient density and the weak wave velocity 


squared. Each of the Hugoniot pressures is multiplied by the expression 

(1 - oe ) and these points are least squares fit for D and S using C as 
a@ known parameter. 

3. Vapor Phase 


The vapor equation of state is a fit to the Thomas-Fermi-Dirac model in 


conjunction with a modified y-law gas equation of state (reference 4). 


The equation appears in PUFF as 


pan(k + (G-H)nt/?] [E-£,(1-expS1- *)3) , 


le 
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where n is 9/06 » H is (y-1) effective for small n's, E. is the 


sublimation energy, and WN is determined from the expression 


c 
nS G(E-o,) 
C, GS, and E are defined above. 


It is easy to show that the vapor and solid equations have a common limit 


of GpE as the limits nel” and n+l’ are taken respectively. 


h, Two-Wave Formulation 


The most current version of PUFF contains an option for a two.-wave 


structure in the solid phase equaticn. See figure 2. 


0 CUSPA ‘e 


Figure 2. Two-Wave Structure 
For u's such that O < yw < CUSPA, the pressure is computed by the normal solid 
equation. For u's such that CUSPA < 4, the pressure is computed using the 


expression 


P = CUSPl + C, (u-CUSPA ) + D, (u=cuspa }* + s, (u-cuspa)? + Go'E 


13 
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Cy > Dy » and 8) are computed in the same manner as C ,D, and S&S after 
translating the origin to the point (CUSPA, CUSP1). 
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SECTION VII 


INPUT VARIABLES 


This section gives a workable definition of all input variables used in 
both PUFF and P PUFF. Those followed with a number 1 are used only in PUFF. 
Those followed with a number 2 are used only in P PUFF. All other variables 


are common to both codes. 


AA(M,1) (i) Constants used in computing absorption 
coefficients. The dimensions are 
om” /em. See section on Generator sub- 


routine. 


ANGLE (1) Angle in degrees of the X-ray path measured 
from the normal of the material. A zere 


value indicates perpendicular alignment. 


B (M,I) (1) Dimensionless constants used in computing 
absorption coefficients. See section on 


Generator subroutine. 


CKP A problem stop parameter. The problem 
will terminate when the maximum pressure 


reaches this distance in centimeters, 


CUSPA(M} Woo -l, corresponding to an inflection 


point in the Hugoniot data. 


CUSPC(M) Equation-of-state constants used in the 
CUSPD(M) two-wave solid phase equation The units 
CUSPS(M) F A 
are dynes/cm . See section on Equation 
of State. 

CUSPG (M) Grtineisen parameter for the two-wave 
solid phase equation. It is dimension-~ 
less. 

CUSFiCM) Value of the pressure in dynes fem’ cor- 


responding to an inflection point in the 


Hugoniot data. 


1s 
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DISCPT Description of the problem. Up to 80 
permissible Hollerith characters may 
be used. 

DA Size of the first zone in centimeters. 

EE(L) {1} Total energy in cal/em® of each black 


body. For problems where ANGLE # 0, 
the total energy for each black body 
should be multiplied by cos (ANGLE). 


EDGE(M,! (1) Values of h*v where the absorption 


coefficients are discontinuous. 


EI{T) ti) Energy in cal/em* corresponding to 
each interval of h-v in an arbitrary 
spectrum problem. For problems where 
ANGLE # 0, each value of EI should 
be multiplied by cos (ANGLE). EI is 
computed in a black body problem. 


EQSTC(M) Equation-of-state constants used in 

EQSTD(M) the solid phase equation The units 

EQSTS(M) > : 
are dynes/em . See section on Equa- 
tion of State. 

EQSTE(M) Sublimation energy of material M in 
ergs/gm. 

EQSTG(M) Griineisen parameter. It is dimension- 
less. 

EQSTH(M) (y-1) effective for low valves of 
p/P. It is dimensionless. 

JBND(M) Zone number of the last zone in « 


material (M}. The JBND of the last 
material shouid be initializei as zero 
for the code, i.e., JBND (NMTRLS)#0. 

In a plate slap, if the flyer and target 
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JICYCS 


JEDIT(I) 


JFIN 


JRZL 


JZPUL 


LOZHIZ 


MATL(M) 


NBB 


NHNU 


(1) 


are thé sume material, the 
JBND(1) = JFIN2. 


A problem stop parameter. The problem 
will terminate after running this many 


cycles. 


Zone number where Lagrangian edits are 


desired. 


zone number of the last zone in the 


problem. 


Arbitrary number of zones desired in 
the collapse region of rezone. Should 
be 20 to 50 for one-pulse problems, 40 
to 100 for two-pulse problems, and 20 
to 100 for plate slap problems. 


Zoning number used in rezone for di- 
viding in front of the pulse. A value 
of 40 is good for all problems. ZPUL 
is set equal to JZPUL. 


Flag used to allow use of rezone ina 
two-pulse problem. Set to zero for 
one-pulse problems. Set to 1 if the 
right pulse is predicted to be dominant 


in two-pulse problems. 
Material name or description. 


Number of black bodies used. Should be 


1 for an arbitrary spectrum problem 


Number of matched values of hiv 
intervais and energy used in the arbi- 
trary spectrum input. If zero, used as 
a flag to select black body energy 


computation. 


toe 
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NPRIN 


NRZ 


NTAPE 


NTEDT 


NTEST 


NZ(I) 


PMIN(M) 


Number of Lagrangian edits selected. 


Number of materials. Ina plate slap, 
if the flyer and target are of the 
same material, then NMTRLS = 2. 


Number of edges (discontinuities in 
absorption coefficients) for each 
material. 


Controls the frequency of printed edits. 
Edits will be printed at cycles which 


are integer multiples of this number. 


Number of zoning ratios used in initial 
zoning. 


Controls frequency of data dumps on 
binary tape. Information will be stored 
at cycles which are integer mltiples 


of this number. 
Number of time edits selected. 


A read check variable. Should always be 
set to 30. If there are too few or too 


many data cards, the problem will stop. 


Zoning number which changes the zoning 
ratio (RZ). The code requires this 
number to be 1 less than the zone where 
the ratio is wished to be changed. 


Minimum pressure in dynes /em* allowed in 
the mesh for material M. Approximates 


the dynamic tensile strength. 
Ambient density in gm/om> of the 
material M. 


Zoning ratio used in geometrically 


zoning a problem. 
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TEDIT(I) 


TIME 


TS 


UFACE 


UFIN 


UZERO 


(1) 


(2) 


(2) 
(2) 


X-ray shine time in seconds. 


Temperature in kev of each black body. 
May be read as zero in an arbitrary 


spectrum problem. 


Storage for hv/KT values for black 
body problems or hv values for 


arbitrary spectrum problems. 


Problem times in seconds where edits 


are desired. 


First time step in seconds. This 
variable is later used for total 


ot? seconds is a 


problem time. i 
reasonable figure for most problems. 
A problem stop parameter. The 
problem will terminate when TIME 


reaches this time in seconds. 


A constant used to smear the velocity 
discontinuity between the flyer and 
the target over five zones. See 


section on Generator. 
last zone of the flyer plate. 


Velocity of the flyer plate. 
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APPENDIX I 


FORTRAN LISTING OF PUFF 


a 


rea ers 


PROGRAM FUFF 
THIS I& PUFF FOR DEPOSITION UNLY FOR THLE AFWL 1604 Rene 


COMMON £A6220)0AC(100109+6) 6B16020) eCUSPI] (6) eCUSPAI6) «CUSPC(6)6 
1CUSPO(6) eC USPG(6) eCUSPH 16) eDISCPTI12) eEE (10) eFOGE (6020) cE LCI De1 OD) 
2ceEOSTC(6)sEOSTOI6) eEQETE(6) sEGSTGI6) sEQSTH(IG) cEQSTN(6) sEGSTS(Od6 
BIBND(6)6JEDITC10) eJORGOIO) MATL (6)-+NOE (6) «NZ020) ePMIN( 6) ¢RHOC6) « 
GRZ620) eSS(BOI VATED) ce TELCIO9O)  eTEVIT( 25) «x( B01) 


COMMON CrP eDTNHeOTNe UCR ee UCYCS 0 UF IMs UPMAX 5 JUOMAMT ¢e UPL ee USTAR ep UKKXG ee 
TPUL. eLINE tL CZHI 2 eNe NES eNJEDI Te NMTRL Ge NPRINONRZONTAPE ONTEOT ePMAXe 


2°OUR 6 TIME + TS e WTAPE e 2OUL 
DIMENSION PCL OCD eCSCL VEC LP eVEL eZMOELDeUCT) 


EQUIVALENCE (€ACeP)  (40( B02) 6G) eo fLAC(1603) 6CS)e (AC(2404) of pe (ACK 3205 
1d eV eC AC (4606) 02M) 6 (4014807) 0) 


ZEPOES COMMON 
oC 1 K21+9954 
BACK) zl e 

CALL GENRAT 


INTTIGLIZE COUNTERS AND CONSTANTS 
LINE SC : 
FLAG=O6 
CC=1e8 
Cl=025 
TWOCL=2,#Cl 
CosgeC7#CoO 
FORC.9=4 © *CCSO 
Tiel 
N=l 
SSC<2065 
OTINS TIVE 
OTNH=ETIME 
HYDRO STERTS HERE 

TIME coop 


SMeMsOe 

PMAX=O06 

Mel 

tLe! 

DETERMINE THe LET SOUNDARY CONDITIONS 
UCLIZUCTIM“DOTINGN(GC294F12) 972602) 

KCL) = 01240 TNHAUT LD) 


MAIN LCOP FOR HYDRY CALCULATION 


22 


Jt rN=1 
Lt Peg IA, 


J*i) 


| TFC GH Pe N ON 1D Ge eG 
LU-LL +] 
q 
q 


56C¢ 

S577 

« se té Ceo TAR) DVebset 580 
Lot OVEMINIF CTI SCHL TN Uae YS Oe?) 590 

TECVC G+) FePet 600 

* fo tbe4/VluF.) 617 
ve ios ec 

1 Fe 630 

CALL ta Ter Cuse yall eh Cue. yell) 640 

‘ eOot DSU Ost be TN® (BP CUF2 FUle te d-Bl GF lest LIS O2Z 6 CSF lve 2iituer yy 606 
Ire. Chet b dv iwber =F) ToekOe)Il 405 

: VOID Ee! ers 
tt MOS* Lr w Ce dei Thee el dD ROO 
OME) ye) Cn 

VOSJ*# Ld =O KCS+) HK ODI ZOD) 700 
VEVGZIV OC 54l}+VOLDIZE of 717 
evVsuTNeeOUs I Ue] ) Te 

TF (DuU4+ted Leettel 3 730 

ee, MOVFH LIT OCMUH CIS GHRCLHCE CUS] DP EMUSVAVG 740 
IF (uC U+# bebe) LaelGelG 7&0 

i* VoIetyere THO 
SUF @ 770 

"’ Cy} ele 7ac 

ia IF CTIMS eae scuw) lel Felo 790 
eas (Ut yee CoD VOSS (yet Pe TNE 200 
Sc 7S 1Pf 41d 

16 TP OSSCK) Lees Felt 820 
i? ECSHPVET CVF L Pee CHD VCS OUR OTIS ai TN) a3c 
a tesle : “ R40 

pa Nv welesvut+l) RES 
CL ist (4d) ABO 

eSRE GT Cl Ove ve whbhiwePoe) 270 

~bFL CoH] VeRO lC IF] PREV PRO 

ALE LGST CkhlelL No wei aM) 890 

CUFL dS (CGF) IK (OS4Ol Ge} Ve Ql ues re eOLLIROVs Ce ee (PL MP2) 7/3 (441)) 900 

CIF) PSPO+ lh ClaHOl ye (CE CUslpyet SLs OLSs7TV 919 
WITOze T1#ONEW 220 

TF (DV) 2 altel 930 

1% DF UDGECN: WA OL TA 940 
5¢ TO 1 _ ' 9&9 

eee JFUCGS0 NEWELL Te 960 
roan | FALL EOST CE CUF1L) eDE UG ePFUTE OM) 970 
LEFURHO= (PF UCGHEP( U4] IS CDFUDGEDNE wv) 9HO 

f€ (DPDRKHO) F2022427 age 

“9 “OFF DS « 1o00¢c 
2d TO Pa ; 1010 

ae SPELL DES ARTE COPORHO) 10<co 
Oe &* QML ae CCPOPED+TWOCL ets (Cue) VeFARESN Gaus OKO i pew E ID) 1924 
lL tvebys PEED 19049 

10.F 

t= CFC +1 PRPMAKD C60eE60c" 106¢ 

wee PVaKMeeC C41) 1676 
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ab WOK ES 


Gao um & 


3 
3 


42 


JRPMA Ks 54} ne 
TF C5 BMH}SK OM] 2? e2 74K ee 
SkamM25K eM] L3cs 
Ma@ubu PPEee 
TF (ue vet) ) 33429.33 Lass 
IF (Net) 30430.31 140 
JPMGX)sPMAxe Ll iie?t 
FFCLOPZHEZ) 3Bise3&eGl ee pons) 
TF o¢g—M5 7423 33.324 30 Peo 
ASTADs J Ly7s 
GO TO 34 218O 
CONT TINUE lieved 
VAIN LOOP ENOS HERE 3.90 

Lear 

USTARSe VE TMM) Pete 


EXIT AND EDIT CONTROL 


z 
an 2 


TF CTEMEH}$TS} B5e27427 
TF (Nw gCVvCS) 36037437 
IF €XCGPMAKj=CKP) Gied?sa7 


ee 
ofa fa 7. fe 


Qos 


DINZOTN4 


c 

oS 

c 

FLAG=L.e 2omed 
WTAPERalLe tree 
CALL EDIT bebsG 
END FILE @& P3226 
END FILE 6 P228o 
RiwinD 4 [sat 
REwIND 6 Meet on 
IF CFLAG) 39.40.39 LAC 
Rt wine 45 Paks 
Stor renal 
CALL DOTF bard 
PAUSE One) 
IF €SENSE SWITCH 2) 42444 pale 
PRINT 439N yeaa 
FORMAT (31H SENSE SWITCH 2 $2 ON AY Owls Pied 14639 
oc FO 328 beet 
LE OXMODE (Ne NTAPE)) 66145.46 best 
WTAPE SL, Lae 
Call ENTT ke 7S 
SO TO 4&4 Laers 
TF UXMOTE (Ne NPOLS.) ) 65447449 14aace 
WARES ¢ OGC 
Chieie BA it ee 
CYCLE ACVANCE Regs 

SROMEMINTE Ce GZEK OM 01 ss 2RlOTNA) 1Sa° 
IF ENJES IT) 69981946 : 
WRITS TADE 4 eNVEDT Tame TIME a 
BDO £3 feleNucCllT ae 
Je2JbolReL) ay SE 
aPiTe TAPE @eJORGC 1) eur CITC sy eFC ens] fae aS 
IF (SDURRTIME) S2eS2. 827 ae, 
Se OMEMINIF Ce Ob Boon se SKE) S 


fog 


14 


ban.) 
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[Ik CTECITCOTII) [865814546 
IF (NTEOT) 55057455 

TF CTIME4O 7 Nee Te CIT¢il)) 56+58e56 
CTNHSTECITC ITIL) -TIME 

NTE OTs: 

So TO =A 

wTAPEzle 

CALL Et 

Tlette. 

NTEOT2S} 

TIMGsT] Me eUuTie 
CINZO TNS MTNGS 

N=N4+} 

TF (OTN) Se EGO 
PAYUSESS 

EO TO ? 

ENL 


165C¢ 
1666 
16706 
+6SC 
1696 
1700 
171C 
1720 
173¢ 
1740 
175C¢ 
1760 
1770 
1780" 
179¢ 
1900 
Hic 


~ 


wo fe 


Ou Fa 


oe 


GENERATOR SUBROUTINE FOR THE PUFF HYDRO CODE 


SUBROUTINE GENRAT 

SOMMON AAG 920) 0 ACL 10010506) (6020) eCUlP1 (6) eCUSPAI6) sCUSPCIUD. 
1CUSPD (6) eCUSPG (6) eCUSPS(6) eOISCPT(12)0EE (10) LOGE (6020) cE LE10010") 
2eEQSTC (6) eEGSTO(6) eEGSTE (6) eEQSTGI6) eLOSTH(6) eEGSTNIG) sEQSTS(6) ¢ 
JYBNC (6) eJEDITOCIC) eYORG(610) eMATL (6) eNCE (6) eNZ(620) ePMIN( 6) eRHOE) 
GRAZ620) 6586801) eTO10) se TEL OC1LN9) e TECIT( 25) «X(901) 


LCYMON CKPeUTNHeOTNe UCR ep UC YCS UF ING JPY AK se JUOMAK! eURZL pe USTARSUXXD eT 


IPUL eL INE eLOZHIZeNeNBBeNJEDIT«ANMTRLS eNPRINeENRZONTAPE WNTEUT pPMAX, 
ESDUR e TIME oe Toe WTLPE « 7PUL 


DIMENSTON POL eOC1 eCSCL PCE CLeVOEED Ce ZME12 eUC1) 


EQUIVALENCE CAC eP) 0 (AC(BO2]) 0G) eo (AC(1602) CS) 0 (AL( 2404) oti de (ATL S7NF 
DreViel AC (4006) 02M) 6 (AC( 4807) eV) 


FORMAT£10ARB) 

FORMAT (8110) 

FORMAT (BELCe3) 

FORMAT (2E1567) 

FORMAT(CIHOTHETITOT =E15e7) 

FORMAT (CIHCBAHLOZrIZ 2119) 

FORMAT (45H ERROR--J]NPUT CARDS AFE NOT PROPERLY PREPARED) 

FORMAT (25HI1 INPUT PARAMETERS FOR = =/1HO10AB) 

FORMAT (G6H ####* THIS PROBLEM WAS RUN ON PUFFSIEOG HK HHH) 

FORMAT (13h TABLE VALUES SOCIOL LGe3e7)) 

FORVYAT (1H06Xe BHNDBBX eQGHNRZC4X oe GHNM TRE 356% e GHURZL EX e SHIZPULS& e SHNPS. 
LINEX eSHANGLE/61100F 1003) 

FORMAT (//17H 27ONING CONSTANTS S/) : 
FORMAT (3 e6H 24:110E1003+8H TO ZONE!I6) 

FORMAT (//45H BLACK BODY TEMPERATURE AND ASSCTIATEOO ENERGY) 

FORMAT (12H TEMPERATURE «5X%e8H ENERGY ) 

FORMAT (1H E1L0e3e6xXeE1003) 

FORMAT (1HOSXeGHIF INSX se SHICYCSSX e SHNTES TSX eSHNTAPE 7x e SHICK PER se CMT SE 
1X eGHT IMEGX + GHSOURs4G110046F£100e3) 

FORMAT (25HOMATERIAL PROPERTIES FOR ABeSX e4HRHOSELOe2eDKe THFROA UH 
11S¢2xKeSHTO Jets) 

FORMAT (109K eSHEQSTCIOX se SHEGSTOIGX se SHLOSTEILOK se SHEOSTGIOX eSHEGSTHI 
1 e SHEGSTSICK eSHEQSTNII Xe GHPMIN/SBE1S5eS) 

FORMAT (1H99xX e SHCUSPI10X e SHCUSPALOX e SHCUSPCIONW se SUCUSPD1Ox eSHCUTOG 
L1OxeSHCUSPS/5OE1S eS) 

FORMAT (///5H NOE&1S5e23u 43H ANelGxXelh 1905 0Ke aM EDGES) 

FORMAT (E37 eS eE C0 e Se E 2505) 

FORMAT (1L41e10ABe/5HO J 6K e2HDOKI3K se LHXILIK eOHERGS1OKe BHCALIOKe 
17HSUM CAL 7X e THERGS/GM7X eGHVELOCITYSX+¢GHZONE MASSGXe lL HUIS) 

FORMAT (14 [3080614 e5e14) 

FQOPMAT (1H1) 

FORMAT (IHCIOXeI GHTHE TEDITS ARE/SIOEINs3/) 

FORMAT (UHCLOMeIGHTHE YEDITS ARESIOIIO) 

FORMAT (21H MATERIAL THICKNESS 2£10¢63) 


READ DATA 
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4a. 


BEG FenmNUeNTECTANUJETNITeLOZ IZ 

TF OCNENG) 29¢35e2S5 

POae Sel Tole byendthelebeleNANy) 

30 TO 32 

FEAL Jet TChryelelel tg) 

BEAL Le tribsCPT CL rel=10e10) 

TF CONTEST) WBeBte%e 

GIAGC Se TES TT OCLD ele) eNTEDT) 

IFONIJLUTT) 38435036 

REAS Se UEOITO) eleleNJEOIT) 

Beau 2 eNuote NP ZC eNMTOLS eURZL ae JZPUL eNPRINe NTADE 
ZPUR FIZPUL 

GES 3eANGLE 

NMTSNM TALS a1 

Emi Ca (NOE CM) eM2leNVTELSE ) 

MEAD .eS57URACTIL) ef (CL) eed oe NS) 

IF OCNMT) 37637636 

ThGG Se ( USBNOOM) MRL eNVT) 

MEAD 2euFINe CNZ(L) eL2leN22C) 

REAL 2edDKe TIME s (QRZ) eMBLeNQZC) 

oxKEDXsR 21) 

CO 34 “sleNVYTRLS 

SEAL 1eVATL OM) 

READ FeRHO() eEQSTCOM) eEQSTDIM) eEQSTE(M) eEQSTG(IM) eEGSTHEIM) cEQSTSEM 
1). ePMING™M) 

SEAL 3eCUSP1(M) eCUSEL IM) eCUSPC(M) eCUSPD(M) eCUSPGIM) eCUSPS(M) 
LOSTNOM) BEQSTCOADSEQSTG(MIZ(CEQSTE (CM) #RHOUM) ) 
NOLUSNCE OM) 

EAC Be (AAC TL ealMel) LOGE (Mele l=} eNCeu) 
wEAD 3eCKPe TS 

GOAC 2SesCVCs eNTEST 

TF CONTE ETARD) 390460039 

PRINT 7 

[TCP 

JCK sO 

SRKOZS 

NRZzeS. 

ANGLE SC OCF CANGLE/35702957795) 

IFCNJECIT) 4344 4} 

SO G2 teleNJEDIT 

YORPGI(IysyEOITCT) 


CALCULATE ALSORPTION COUFFICIENTS 


TF (NANG) 44¢50 e446 

CO 48 VeleNMTRLES 

La I | 

30 48 TeleNHNU as 
IF (EDGE (MeK)=TILIT)) 46046047 
<2K+)] 

6C TO «48 

Al Lele) @=RHC(M) FAA (MEK) #(TBL OT) HR (MOK) )ZANGLE 
CONTINVE 

ELTOT8 © 

30 49 TT teleNRNy 
ELTOTSESTOTSE Pde TT) 

$0 TC EM 
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e 


ft. 4 


wife 


2680 
L600 
z=700 
710 
2769 
2735 
= 74D 
27295 
«< 760 
«776 
“790 
2790 
€800 
2sl¢ 
2725 
2130 
£2446 
ees 
2160 
aero 
2kB0 
PSO 
2900 
eOlO 
e206 
A439 
EIGd 


aes 


EITOT=0, 

00 S&@ M2z}eNMTRLS 

0O S54 Ls!eNSB 

Kel 

00 S4& 121.109 

TF (CE0GECMeKIKTBL (1 *#T(L 9) 52452452 
Kak+} 

GO TO Si 


AC(L eT eM) ZeRHOCMIFAACM =< eC TBL CLI #T (Le 8B OMe K) SANGLE 


CONT INVE 

00 S7 t=z1eNe8 

DO S® Im1e¢189 

TF ¢(1=-$9) 55255+56 
El(Lel2EE(us#e01 


GO TO 57 
FIldine 1 SEE(L 24.001 
CONTINUE 

" COMPUTE Ox USING ZONEING CONSTANTS 
LZ} 


OO 61 JeleJFIN 

TF (CJT eNZ(LZ)) 60489065 

LZ2LZ+1 

OxsOxK#RZ(LZ) 

xCJ+1)2«09)4Dx ‘ 


ZONE CESOSITION 


Vol 

CO 71 yetesFIn 

ESuUMtSe 

IF (Je JANC(M)) 6346276. 

M=V41 

fF (NHNU) 64+67+64 

CC 66 [el«eNHNu 

IE CEM pel elel-20) 66265065 

EL7FeE Cpe lL VBC LemEXPF COC CL el eM HI XC US] eK EID) 
Elchbetlsyellitel rElz 

ESUMBESUMFE TZ 

CONT INVE 

SO TO 7% 

CC 6° U21+eNBB 

OC 69 T2tel tS 

IF CEt(uaelbeleE"20) 69+68+68 

CpPaT (Lele CL ermEWMPOF (CAC (Le leMPH(K( el yew)? 
El(ueleli(lLelreEl2 

ES UMEESUM+E TZ 

CONT INVE 

G2 yel  eEOUMHG, LBBEFZREO(M)/S(X(J41 2 XK6 991 7S0UR 
1- (SSC ys+beleLlAsReOty) FE. 7Le 7) 

CANTINJE 


CLEAR STORAGE TO ZERO FOR HYORO 


Co 75 te#1+e569r 


Alc ele 
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2960 
2g7O 
aes 
2930 
3906 
36:6 
3020 
3030 
3040 
305¢ 
3055 
3070 
3065 
3099 
a1ce¢ 
3his 
3125 
3135 
3140 
3150 
3166 
3170 
31A6 
3190 
3200 
3210 
3226 
3230 
3246 
3250 
3260 
3270 
3286 
3290 
3300 
3310 
232C¢ 
3336 
3340 
335¢ 
3360 
3375 
338¢ 
3396 
3400 
2416 
342C 
3425 
3aae 
3455 
3465 
247% 
ine 
3496 
350° 
35130 
352% 


Best AV giiapie COP? 


yun 


mae 4 
74 


SET YUE ZONING FOR HYORO 


Me} 

990 FR seeesF tn 
VOUitleasimoiM) 

PMCS SCKCId=K CURL PIZVEOI) 
TF (5-URNDOM)) 75074075 


Merve] 
CONTINGE 
START INPUT EDIT 
PINT + 
PRINT Te(BLECPT(K) em elel 
PENT 4 


POINT Te (THLE Tel a1e109}3 
PRINT Ll eNSSeNRZOeNV TELS oe URZL ee IZPUL eNPRING ANGLE 
ERINT &.81TOT 

OR INT aL IFHIZ 

OnINT 12 

PRINT 13e(QZ0T reNZCT eT 2b eNRZC) 

POINT 14 

POINT 156 

PEINT 1506 TCP CEE CT De LE1eNHS) 

POINT 1 PeJFINeUCYCS eNTEST eNTAPE eCKPe TSe TIME ef DUR 
PE CNITIITY 770770 76 

QIN Pret JEOUT OD elslemyeorty 

PE CNT UST) Tee Te Te 

PQINT 2HeCTEOIT OL el eleNTEDT) 

Jere 

RINT PS 

27 BG VeleNVTALE 

TF CUNT OM)) 21e3C081 

Jslegrds 

GO TO 27 

J OC BUANTOM) 

TRENGSx 592 =k C81) 

SCOINT (, OeMATL OV) eRHO(M) CURE OURS 

PELE Ine 

PUINT 2PQe THENE 

E-INT 19¢fQSTCOHr) eESGSTOCM) eEGSTE CM) eEGSTG(M) sEQSTHIM) sEGSTS(™) eb OS 
TN OM) Ge OMEN ON) 

PE OO SOS OMY Y RGA eA? 

RINT 24 6CUSP1 OM) oe CUSRA EM) eCUEPCIO™M) eCUCOD(M) eCUSPGIM) » CUSPS (14) 
NUE OSNOE (4) 

POINT PLeNCED 

PEINT 220 ALI Meld eS (Mel) eEDOGE (Mel) ef B1eNOED) 


CONT INGE 
ITART OEPSSiTicNn EOit 


POINT 2PAeCOISCPT Ct de betel O) 
M=] 

SUMCAL ECs 

CO 92 YH2eIFIN 
Boss 37 (9) *3S0UR 
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7 CGH OE NOT OM 4b Be AT 
Ascte] 

pe {73 eR TET OM Rd eo 7) FIST g rte THO 
TE CGETLO)Y BQg fg GN 


Jen Jel) 
sFOGHeRHD (vy #Ox 
“ T54PAe1L eH 7/4, 186 


LeTuMCAL OAL PA 


Fix 5d) Gel es 
EQINT MMe CT PTCOTE Pye fete l™) 


. ie or 
Bo, TAL Pte te Dt e me 
JET Is ETN 
sryt ote 
ord nal an Sd 
cass 
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PH tema KOS) OE EG™ pe CAL Phe  UNCAL ght 


N 


eUCIIMEZM CI ou 


SPL cet Tr & Or tee Mure Hw e Cols eons 


Gar eS EE See eon 
aot 
AMD AACS OECD AAC ELGG Feed et Oe LO) eCUSPI (6) sCUSPAIG) eCUSPC lO). aren a 
TE Cy eT UEPE TE De TOT MTA DA TIESCETOCI2S SELECT OVAFCGF (Ge lO AE POLNIG +) gen 
See TST CE De GSTT EB) AE USTE (ED AE OSTC(6) HEGSTHIB) CE GSIN(IG) eCOSTS(O)e 4.140 
Ba CA pe CLT OD  eITE AECL OV eA MATL OG) eNSE CED ONZO2%) CEMINGG) ePHS (Ode 4220 
iO, 7 Ye Co vel OEM eT LOD Tre TEODOR eKCAST) 464 
aave 
Se a teh eg Tite Teel eur. aul INe RP MAX e UPMAKT e UR ZL oe ISTAR OG UXKOe .’ “4 tek 
Pon et the aL LOM Yate’ tee OTT eM TELS eNPRINeENRZONTAPE ONTEDT ePhiAxke 4.390 
ea ge The TAO gg TOL asco 
4ald 
. TOM SOD De Oda TSOL de Cl deve FMC LeU) wee D 
aa 
oe PV LEME CAC oe CASON DAG oe (ACCIGIS elare (AlC(CeGOG) eb CALCOT. 44a", 
ee OOO EA TB a ae LAC (4RDT) ou) ic Marc, 
ee 4685 
wa hl 
20 at oF [Pe See lI) 44r5 
aacct 
REPOS lee AD CF MAIN PRESSURE PULSE @5n9 
' asic 
ClANCH2 6 asSer 
vTapcety, 4330 
TE (GR PAK RUPE LL MSR ZL MUP AXE) B7eE7 ee G£4aCc 
uy} GEBS 
“Se 6560 
Cf & J2leuySTAP 4570 
YE (UR JANDOM)) G@e3eu asae 
Ma. +i 45SC 
To CRHT OMI EV C4 lH be 2) Fe Sek 460C 
JIVeN 461 
CONTINUE 4620 
TF C2PUL) 2B e4Re% 4639 
wv = JP tAwWes 464C 
wD Awsw ( JOMAY } w6&ot 
crevPus : GE69 
T OVPL SE MYPL4UCUY te SH C7 (CUMS +2 Cu) ) GA7C 
TE CUCUMeE)) PFet Pere 46> 
© CUMet? Plelsell 4609 
suey el “790 
we Te. a?ic 
vite JOVEK 4G “w7cd 
Syms. 4730 
RYE SS My DHE IAI SE aM ETM CIM DO ZM (MOTD) 4746 
re 4 JMR TAN) 14eteeiFf a7se 
a] o7é% 
€.: rs 4 PTD 
DT Ye lek he Yer we Ped 
YE J OMAK 745 
Pek TE RIN CEVA) ZOE aoe 
ge yottayr “- 19 
wrel aS 
fr Or ed a es ee oh dod 4-7t 
3 aeud 


Pe GL pew seb yaew ed be edited dS 


1e@ 
is 


eo 


28 


29 


TF € J@ JPMAKM<220} 16446168 

tPF €4-5xKO) 48248.20 

Cale FOUT 

SCAQNCHE!L 

JF INGs UF IN 

IuMKGsS 

DO 22 MOKal eNMTOQLS 

TF €J<t-JBNOCMCK)) 23+21,522 

JCK eMCK 

CONT TINUE 

MCK ENMTOLS 

ME=MCK 

TF CCK), 25.25626 

RZRzC sSHCMOCUAL PERC UDI RK CUR 2) SEK DEKE eH 27) 
FLee tum] p4R ZRH CHE Cet) 

WLW O-1 PER ZROCV OS) Visel)) 

SSL=SS (Ue 1)+RZRHCSS CUES C11) 

IF (U=JGND(ME}) 264274526 

RZRECKC Ye) pame BECK CUP eK Ce LEIS OKC CHE] eK ER] )) 
ER=E(CU+L4#+RZROCECUIMKE (G41) ) 
VREVOU4})  4QZRe(C VE Jde VO Uel)) 
SSReESSCUF1LI+RZRe¥(SS( USSF V41)) 

IF (4CK) 29429428 
RZR=afeSREKC IIE KC UHL PL OX OCIA K ER 20) 
ER=E (CU eRZReCE (Cu )-E( us) )3 
VREVISPERZREIV OS) HV (eb) 

SSR2SE (YF RZRe¥CSS CG) H$SSt uel 

GO TO 2g 
RZR=CeSHOKL YIN K (Gel pietiax toed dextuet}} 
ELSE (Jy 4RZRO(E (CU) MF (oe) 
VLEVEJIERZREI VES evi seh?) 
SSLBSS(JIFRZRE(SS( UP MSS (G4) 3) 
RROX(U) 

Mise Sei xt uryeXtueld? 

ZMReHixRexKL} /VR 

ZMU at KReXL I SVE 

VREUCS) 


Us CC7MC JPW ZML RUC Ue Pe ZM EK 2MR) HOLE ZML+2MR } 


DENELSe SVL 

CALL EQST CEL eDENe PL ame} 
DENF le /VR 

CALL FGST CER aDEN,yPReME) 
W=GtI) 


OG 31) NME BUDO eUFIN 

ITF (NME @#JBNCOCS1)) 3143031 
JAND (YT) =NME +} 

Jl=Jl+] 

CONTINUE 

Ixus 

Ja uF IN 

K(Uel)aexi) 

Utd+ldeutu 
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ws ow ld 
aon 


aR 
49 


Pe ee we eee ee 
Vesta) 
ae ee ae ee ae 
BPlyur1lysPy) 
WC4+1702400C) 
LG6se} 2eSSty) 
(lo be° 709) 
Jzunrl : i a ae eee 
[Fo tyeux™K) 33¢33e32 
xf Je] )=xKR 
KOSS XL 
UCu+ dy sue 
UOCSPSUL 
PVCS) YSZMQ 
CM CVS 
VtutFlIsVve 
YVOSSVL 
pa eo Oe 8 7 
ASP EFE 
Hy y+! >} pron 
PCI) FPL 
VOU+>D)=29R 
VOIVSSL 
CoC us pete 
"ee gp]C Se 
LC J*b) SSR 
SS CSP ESSE 
TE CNJEDITT) 47 e37e34 
CC 46 TPS leNJEDrT 
IF CSPVEDITOLITD) 35035436 
YECITCIPISYVEDITCLL 41 
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PROGRAM P PUFF 

THIS 15 PUFF FOR PLATE SLAPS ONLY FOF THE AWFL 1604 ## 8H 

COMMON €S$(5900) «CUSP1 (6) seCUSPAI6) »-CUSPC (6) «CUSPD(6) + CUSPG(6). 
1CUSPS5(6: ~ITSCPT( $2) +6E(500) eEGSTC(6) eEOSTDO(6) sEQSTE(S) sEGSTGI 
26) sEQSTH( EG) sEQSTNIG) eEGSTS(6) eYBND(6) sYEDIT(C 10) eUOCRSI19) eMATLI5) efi 
3Z(620) 4P 1500) ePMIN( 6?) Gt S00) eRHO(6) RICZO) ee TEOI TEAS eUlC TOG eveer*C)e 
4x(500)«Z2Mt(500) 


COMMON CKP«DTNeDTINHe UCK ep UC YCSe UF INS UF IN2s JPMAX se IROL pe US TARP IKKC eT 
LUL eL INE eNeNJEDIT eNMTRLS eNPRIN «NR Z eNTAPE eNTEDT cPOMAM eT IVE TS eZTERC. 
2wTAPE.ZPUL 


ZERCES COMMON 


OO ft K21.4279 
CS(K)tCe 
CALL GENRAT 


INITIALIZE COUNTERS AND CONSTANTS 


LINE=C 

FLAG20 2 

Co=1es 

Cle0e25 
TwOCle2,#Cl 
CoSGeCO#co 
FORCSO=4e*COSC 
biel 

Nel 

OTNsTIME 
OTNH2T IME 
JSTAR#JANO(1)+4+3 


HYDRO “TARTS HERE 
TIME LOOP 


SK2M200 

PMAX20. 

Mz} 

tl=1 

DETERMINE Tre LEFT SOUNDARY CONDITIONS 
UCLIZUCI-OTN®(Q62)4P(2))72M(2) 
X¢D)=XC1L)4¢0T NH RUC 1) 


MAIN LOOP FOR HYDPC CALCULATION 


JFINM1 UF IN@1 
DO 22 JeleJFIN! 
QOL D=zCtU+)) 
POLD=Pty+1) 
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UES] PUCIF+1LHDING(P (9427400 S+2d Plot) aSCUtl ISO ZMC 4d POZO 4. 1) 


IF CA0SF (UCI) Pel eERH3) 5¢506 
UCU41)9=000 

KOCUFLI EX CU41 CD TNMAU (Ue) 

DvsU(Je] euUtl) 

VEIL SOXOCIF1 Ie xX (IIS ZM (41) 
VAVGE( VC SFL I FVOLOISE 60 
OVSDTNH#OULZM{ J+] ) 

IF (DUt}e) 7468.48 

SOIFt) SCDURCCSQ-CL#CS (941) ) #DUsSVAVS 
TF (QC U+bdate) BeGeO 

QCI+1) 2000 

Duel eO 

CS(IJ4172008 

DNEW2}e/Viut) ) 

EOLOZE (U4) 

CALL COST (EC I+] ) eONEWeP2 eM) 
ElSEC U4] )-P (U1) #OV 

CALL EQST (E1eONEWePle™) 

EC J+) ) SE CIF) (CO 24P (4) 24005412 +COLD) HOV 2eH (Pl HFC SPL 41)? 
POJ+1)2P24(P2—-P1)#(E 6441) -EOLO)I/SPOLD/OV 
OLTO0=2.CN1#0NEW 

IF (Ov) Jdelleld 

DF UOGZONE W+OL TD 

GO TO le 

OF UDG=ONE W-OL TD 

CALL EQST (ECU41) eOFUNGePFULCS eM) 
DPORHOs (PFUDGHP( U4) DIZ CLEFUDG-ONEW) 
[IF (DPDRHO) 13413014 

SPEEDS, 

GO TO 16 

SPEED2SQRTF (OPORHO) 
SK2M1=(SPEED+TwOCL#CSs (41) -FORCSO*#DU I 4 OKO U4) =X)? 
CStJ+1)=SPEED 


IF (PC U41)—-PMAX) 17616016 
PMAX=P(U+1) 
JPMAXK2 J+] 
IF (SK2M-SKEM1) 16018019 
SK2M25K2M1 
Maul 
IF (UC341)) 22420022 
IFCJMJETAR) 22021021 
JSTARSJ 
GO TO 23 
CONT INUE 
MAIN LOGP ENDS HERE 
JSTAR JF INM 


Exit ANF EDIT CONTROL 


TF ¢TIME]=TS) 24+26:26 
TF (N@uCYCS) 25+£€6026 
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FEO OX CURPMAX)-CKE) 30426026 
FLAGele 

WTAPEAzALe 

Caulk ToT 

END FILE 4 

END FILE 6 

REWIND @ 

REWIND 6 

IF (FLAG) 28+°29+28 

RewIND 45 

Stop 

CALL OOTF 

PAUSE 

IF (SENSE SWITCH 2) 21.33 
PRINT 326N 


FORMST (314 SENSE SwITor 2 Im Of AT CYCLE 


Go TO <7 

LF UXMOOF (Ne NTAPL DD) 39434025 
WTOPE Sle 

C4tL EDIT 

GO TC 37 

IF ¢(XMODF (Ne NPRIN) D 2743637 
wTAPEsve 

CALL EDIT 

CvCLe— AGVANCE 


SKEMEMINIF ( 08/582 61 e2#OTNH) 
IFC NJEDIT) 38240038 

WOITE TAPE GeNJEDITeNe TIME 
00 3Y LlzleNJECIT 

Juss JEOITCI) 

WRITE TAPE GB SYORGI1) eUEVITED) Plt) 
OTN=D TNH 

CTNH=SK2M 

CALL RIZONE 

IE (TEDITCIID) 4546504) 

IF (NTEDT) 42444042 

LF (TIME+OTNR@TODOITOCLI)) 45145045 
DINGS TEDITOLI A TIME 

NTEOTSC 

Go TO 4& 

aeTAPEZle 

CALL DIT 

Pislt+) 

NTEOT=1 

TIMESTIME+OTNH 

DTIN=OTN4¢O TNH 

N=oN*] 

IF COTNH) 4604602 

PSLSEe. 

6c TO 2 

END 
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SENFRATOD TUSBDOUTINE FOr THE PUFF HYCRO CODE 
SUSROUTINE GENRAT 


COMMON €S$(500) «CUSP1 (6) eCUSPALB) se CUSPC(6) eCUEPL (6) eCUSPGI6). 
1CUSPS(6) eDISCPT(12) +E (S00) sECSTCLO) cEGSTOIG) SEGLT (6) eCuUSTSl 
26) VECSTHIG) eEQSTN(IB) eEGSTS16) eJONOCB) PVEDT TOA eIOR ELD oe MATL Ce da” 
3Z629) «P (500) ePMIN(6) 6 OI SOO) sPHOCEeRZEZOC) Ae TEDIET OOS ver OF av W416 
4xtGlCO) «274(590) 

COMMON CKPeCTNeDTNHe UCK AUC VOCS eUF IN UF INS ep JPM AK CUE TL eS TAR gp UKK TE 
PUL sL INE sNeNJECIT aNVTRUS NER Tee ete TAOL eh TEDT RUA ee TIL a To erin te 


ewT4Pe,.ZOuUL 


FOZMAT(C1OA8) 

FORMAT (8110) 

FCRMAT (8E10e2) 

FORMAT (46H ERPOCR--18PUT CAPDS ARE NOT SRORPESLY PPEEAIED) 
FORMAT (2EHLINPUT CARAMETERS FCK = eT HOLCARP) 

FORMAT (G6H #448 THIS PFOELLMY AD RUN ON CHOUFFSIESG teeny 
FORMAT (CP HCSX ee GHNRZC4X eGHNMTRL SEX ee SHRI LEM OHS ZRULS& et SIPRINE Xs 
LGHUF ACES e SHUZEROSX«ASHUFINE/SELIOCZE Loe Fe 110) 

FORMAT (//17H ZONING CONSTANTS/SZ) 

FORMAT (3%e6H RATIOEIC.3+8H TO ZONEI4) 


FORMAT ¢€1HOCSXK ee GHIF INS oe SHICYCS5S Xe SHNTESTS& e SHNTADSE 7X 4 BHOK PRX PHT LE 


VX eGHTIMEsQI110¢3E1C 23) 

FORMAT (C2SHOMATERIGL PROPERTIES FOR AGsSXsGHRHO=E10—30° Ke THFRO Ur 
PISe2eKeSHTO YsIS) 

FOL MAT (IHOOK eSHEGSTCICXK eSHECSTOICK eSHECSTEILOCK se SHEGSTGIOX eSHE USTH! 
10K se SHEQSTSI OX eSHEQSTNIIXseGHPMINSSELS 25) 

FORMAT (1 HOGX se SHCUSPI 10K se SHCUSPALO Ke SHCUSPC IC eSHCUSPD1IOXeSHCULES 
110xXe35HCUSPS/6E1565) 

FORMAT (LHL «eI CAB/SSHO 4S GX eT HKAK se BHVELOCITYSxX*eGHZONE “ASS//) 
FOOMAT (16 13eR3E1Ge%) 

FORMAT C1H)) 

FORMAT CLHOICKsIGHTHE TEDITs ARE/ICEIOeS/) 

FORMAT (CIHC1LOX.IGHTNE YESITS ARLSIOIIO) 

FCRMAT (21H MATERIAL THICKNESS 2F1Ge3) 


PEAC DATA 


REAG 1eC(DISCPT CI) el ete1O) 

READ 2«NRZCeNMTRLE eJRZL ee UTZEUL NPE INeNTAPE FNJEDIT eNTEDT 
IFONTEOT) 21021026 

READ Be ( TEOIT(C I) sli eNTLOT) 

IF CNJEDIT) 230223022 

READ 2e( JEDIT CL eleleNucolt) 
ZPUL = JZPCL 

NMTENMTRES<1 

READ 3+UF ACE «UZEROeUF ING 

JF INQSUF ING 

IF CNMT) 2002544 

READ 2a ( JBNOOM) Vel eNvT) 

READ 2eJFINe CNZOL) eb =] eNRZC) 
READ BePxe TIMES CRZIMdavMELeNQIPC) 
DXSDK/RZO1) 


DO 26 VeleNTALS a 
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REAR Le sATL EM) 


ROEAT BePHOCM) «LUT TOOM) eLQSTEM) el GE TE EMD EQUI GOM] GFR THER gt a FF poo 
DP) eR TEN OM) 1ee 
READ 3eCUSP1I 6M) sCULPAIM) «-CUSPC (OM) sCUSPD OM) eCUSPGEM) eCU PSM) ae SO 
EQSTNIM) SE QSTCOMISZEASTGIMISEC™ QSTE CM) PROM) ) Pa? 3 
REAM BeCKP.T. 2 
READ 2eUCYTSeNTETT es 
TF CNTES THIN) 27 ehSe07 “¢ 
PRINT 4 are 
STCe ee 
CONTINGE ae 
IF ONJZOIT) 32431024 uo 
Cec 3¢ ls=iteNJEOIT 2am 
JORGOCTISIEDITE I) ~ 200 
JCK=9 ee) 
Juxdgen we 
NP? =-5% 125 
“em 

SET PLATE VELOCITY ok OG a23 

-. ‘On 

DO 32 J=leJFIN? CLIe 
UCJPEUZERO chee 
CONTINUE are) 
UCJFIN2Z) SUZFRO#( 1 e+UF ACE) 26 «400 
UCUFIN2Z 41) fUZEROSUFACE c41d 
UCSF IN2 42) 2UCUF IN241 97207 Lact 
L459 

COMPUTE Ox UTING ZONEING CONSTANTS 440 

: Gun 

Lz=! ats 
00 %S ysleJFIN 479 
TF (Se PeNZ7OLZ)) 34¢33034 24FO 
CYL Z7+) 249C 
OM OK#ESZFILT) aloe 
MC U+1) 2x0) 40x aeisS 
eevee 

SET UP ZONING FOR HYDRO ieee hes 

te 

Ma] USS 
Of SF Js2edFIN - 6S 
Velez beSRHO(™) lw 
rte pee Yee KMCUNKlLIIZVE9) are m 
TE (Ge UENDOM)) 37636037 vo4g 
Ys’ eo} . PAOD 
CONTINUE 610 
~ AOS 

START INPUT ENIT ate MY 

coat 

OErINT © .€.0 
ESTINT pe SIPECe TECK er el e149) COED 
PRINT 6 ; P E70 
SHIN T PeNRZC ONUTILI eIPFZL oe VZOUL oe EOIN SUF ACE suUZEROe UF IN? “€30 
Repel A AD 
SENT Se CSP OPINED eT el eNe vr’ 709 
CR INT PTF Lee CVC eNO TeNTAPE Cre Te TIME 710 
IF CN GT OET) “4 4 35,20 OMS 
ere ae bl CPT OD Pe LSLeNJEDIT) L750 
44 

oot’ 

F 1h Gov? 

(noe 
Wake 
ran ‘a 


Pan rae oe ae nanan a IE ES AEE STE NE CRESS a! alae estes mete 


an 
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see pet rare TOOT OP be Psd aNtTl 7) 

Jeist 

OTEST 16 

3C 4? ve leNMTELS 

TP CUSTOM DD G24 ea 

JUOSIEIN ; 

S50 TO 44 

Ue Pe JPEND Ov) 

Tre NGex (U9lD eX CIAL) 

ere PNT Pp eMATL (IM) EHEC) a UP] er 

ORINT 19¢THKNS 

OR INT p PoE QETCOM) ESOT) LST TE OME LT 
DTN OM) eP™uMING YD 

PoCTUE PSL) ) G6 4606% 
Pps. pReTE PLM CET PA ET eT USEC OA 
oe Ps 
CONT DS GE 


“TAPT VELOCITY EDIT 


POINT poe C(RPTSCOTCPyestelbel) 
JSOKSIF ING +3 

OC 4% YsleUsDXx 

POINT 156 Ve KC FP eC IIe ZA) 
IF (XMONF 65450) ) 65944804? 
POAT pae(OISCOTCh re betel) 
TANT TINIE 

PRINT 146 

RE TUR 

Fun 


THOM) OE GE 


SECT CULPSOM) 


moEM be Te 


soe eet eee: moet oe eee & ‘ ye ok 
ere 
Op DQ LTINE CE PONE s07¢ 
apc 
zooeson, £7 S002 eC USK (hyaTU PhS CUEOTIG) eCG PME eCUTIPGO)6 $0Feo 
pr OTE) AMLTECPT OLD) GE (09) 6G TH (OD CEGSTS (6) FOSTE LE) HEQSTC( s1Qd 
DE 96D THUG) CECT TINGE) oH GATE (6) eI CED edt PET CLAY eJORTOVO) MATL OM) a- <217 
ge ori 19 ( SOD) ePMINGG) QU SOT) THOCAD AIP EZOD Ee TELTT EOS D UCT ATV AWE OTD pad Urey 
ace Tye TOS) p20 
BX Cus Le Dem TN a DIN e ICK ICY aE ING SF 14? 9 JOMAX o JATL ed TALC IKK 8 I7E apar 
POL eL IME eNeNJECET AN TELE aOR Ie NEZ TAPE aNTEOTePMAK 6 TIME 6 Tt eure tues 3132 
Da TAPE ea ZR : re are) 
» “370 
FORMAT CFLS e2ettIO? as 
F a ee 
AOPONF Gur As OF LIN TAT SSURF PULSE ese | 
Ke 
SCANT HEV 4 Dake 
WTAPKsle ae 
IF (CUPMAKeJZPUL<YRZL O10? BSeRSec ut 
Jvel os 
Ms] ~ 07 
00 6 yeleJSTAR TS 
IF Cue JINndOM)) 4e2ed artia’ 
MsM+] teat 
Te (RPHOCMP HVE LS 222) Fee 7 4o0 
jVteu 31D 
CoerT (sak aire 
IF CZPULY G6046e% oc 
JvryoMAwet ’ Lah ig 
KUM EK SKC JOMAX) gee 
Ev YyOLE” « 16 
EMYPL EEMVBL SUC US DH eSHOZ CU eb eZ (IID £28 
IF cit utely) Lesblelo a 
IF Code gy Tleitell al 9 
SMe of Beet 
ANTS © “ay 
JM=JEVAK4G Fe 
Syene tan 
EMYDORSEVYPREUOCIT DSA HOP CIA VA ZL CIO] ) ed 
IF CSAS TARY 146616018 Aye: anes 
vaste! , ate 
sO TO ia Bets 
euyePs = 'VPL+EMyrs 140 
AT POSEY PE SEMEX vee) 
DP Mwse ee NTRP eC (JEYVAM I Z79UL eae 
Je JPrar ay toe 
JeJ+} gets 
IF CU+TeuFIN) 17046046 atte 
IF (XC 0) =k CUNT DRRZOX) 1B61Bel9 17.45 
{© ¢ JeyOMaxKeZo 16246046 awe 
TF te KKoO) 46146025 onée 
CALL OTT f279 
SCONCHELe Gee 
JFINOZIE TN. wae 
JuxdeJ 2OOL 
DO 22 MEWS ANETSLS 1BY> 
46 


PO gate ar (6S 
nee Oa 

COM TPM 
Ch ENE TOLLS 

eo eMC 

TE CICK) 250256276 

REQ Co ROME URL DEK UD DRM CURD SOKO SI AKL Hed) 
TLE CIJH1 POP 7ZREOCZ ESD HWF CueMdd) 

Vi eV Sel eR 7ReE CV EC-V Oued) 

I" CUM JENO CME) ) 26027426 

TECK OIe) pm e BHCC eK CIRM D DDI EXC) eK CURD DD) 
“RSC l Use] VeRZNeC EC UINL (U4) ) 

VO=SVO SFL eR 7ReECVEO HVC UST)? 

IF ¢utey 29,29,28 

RIDE (Ce THOKCII RKO INK DDI ZONE IIR KOH)? 

PRE CUP*AMZRHCE CUVMHF CM d) 

VISVE UP +RZQHEV OUI REVO Rl)? 

6c TO 79 

PIMs (ae SelC KCI C NHL VIS OKC Ie] vex CURD) 

ELSE (UpeRZRH(ECIIMKECU41)) 

VL VC UI497Q28( VII) HV EC UO1)) 

KO=K OI) 

XL= eS ROK CUD eK CUM) 

ZMO=(xIe XL I SVP 


ZML = OxR=XLIZVL atta, 


UP sUCd) os a 
ULs€C24° CS 2M FUC Ie 11407052 © 7MP) 80) ISO 2ML +7 MR) 
OFN=le/VL 

CALL EGST CEL LDENSPL EVE) 

OENS1e/VA 
C4LuU FAST 
Gs) 
Greats) 
csueC$(y) 
ese=ceiy) 
JOOrS 
JL=omMe 

CC 31 NMEsU00e VF IN 

IF ONMEeJENOCJIT)) 31430031 
JIND CIT) =NME +} 

JPsJiel 

CONT INVE 

IXKES 

J=JFIN 

KO J*)L =x) 

UCJel sul) 

ZM(IJ41)e2M(09) 

VeJelreveu) 

EC J+ pe) 

PC J+1) =P) 

9941728005) 

CS(J3e)2CStu) 

JeJq=-1 

IF (€J~UxKx) 
KXCUe1)=sxR 
KOI) 2xL 
UCUel) sUR 


CER eDENe EO ME) 


33633.32 


47 


SOU 
467C 
HObD 
by aan] 
wTOS 
2716 
37.0 
374 
3740 
3753 
2760 
AV7S 
ITS 
27aC 
3800 
Aig 
400 
JV30 
24ad 
2RED 
jAS69 
Reese) 
Hag 
a 90 
oOQ 
4910 
2400 
3920 
3960 
3930 
3960 
JIY7TO 
398¢ 
3250 
“occ 
4010 
40<0 
4040 
4040 
“Of 

e16C 
ao7e 
@On€ 
4095 
aioe 
oie 
ree 
aire 
a147 
ace 
416? 
air 


4a41reg 


44 


pet 


- 


VOOdSUL 
de ae }=Z7MR 

26052 27ML 

VOJelpsyn 

VeJ)evE 

Cluel ser 

SO JV2EL 

2( Jel )=PR 

PO) SPL 

QCu+ly=9R 

Qu) =OL 

C5441) =C5R 

ROO ECSL 

IF CRIOVIT) 3PM are 2a 

20 36 TfflseNJEDIT 

TF CJ@SODITOCLED) 35535036 
JLOLTOLPYZSFOITOLL eo! 
CONT DAE 

IF ( JeucTAR) 2B 438,29 
YCTOESsUOTADS] 
SPINS JF ING! 

TF xO Sp =X (G1) =-RZDK) 40.24.24 
TF €I"UxKXKO-20) 61445,45 
Jzs+2 

IT (COME IN) 62.65.45 

TE (GMb JOND OME) ) 42460463 
Jcest 

GO TO 78 

SIoME as 

JoKs] 

SO TO 76 

YJRIZEUFIND@ FIN 

EPIINT Le TIME se eURZ.UCTAR 


REZONE BEHING “OIN DEE ss 


TF (URZL) B5.AS.47 

TE UNMNR 72102) A886 ,an 
JVs] 

“ay et 

YF INOS UFIN 

YF INSS UF ING] 

M=oVin? 

KAPTSe 

J=JPMGx 

Jz} 

TF CPCS SPMAK=62) 504560249 
YPLC ZS 

TF CSJPLCeUVevP2L) 31¢51.58 

re (JOVEK@SELORKI2>uUL ) 2052053 
TF (€37) B5,AS5,R4 

vI JPL +] 

DS SG MMS LaNMTRLS 

TFC JPLC+1-53ND CM) ) 5 055054 
CONT INGE 

RAMS AIA T ILS 


Mann 


48 


70 
71 


74 


7% 


DO 56 MM=zl«NMTPLS 


TF CUPMAK=JRAND IMM) & 7457056 


CONTINUE 

MM=NMTALS 
RZRARZOX# e S*#V ( IPMAXK) HQHO CMM 
KRZs@] ; 
GC TO 61 

KPZ=1 

00 SQ MMe] .NMTBLS 

IF (3V41*JRPND0MM)) 697469689 
CONT I NUE 

MMSMYTELS 

Ms lM 


RZR=(KXCIJPLE =x CIVIIZFLOATF(CUR2ZL) 


JzIVel 

TF (SCRANCH) 624662464 
CALL FOOT 

SCENCHEle 


IF (XC U4+1)=2XCYI-RZR) 60079, 7° 


TF €4941—@JBND(M)} 65079065 
IF C(NJEDIT) 68068266 

90 87 [l=leNJEDIT 

TE CUF+1TeVEDITOCLIID) 67579467 
CONTINUE 
ZMINVELeSC2ZMOC 41 142M (42)) 


VCIF+D I SOV ESF) HZM C541 4V 0542) H2M( 542) DRZMING 
EC J+ bre CE CSF) PHAM OC USL PFE C2 HZME 42) HZYING 


DEN=1le/V0CU+1) 


CALL EGST CE CUFT) CEN PC Ueda eM) 


UMS CZM(SIHUCIIFZM (CIS) DH CUL ID CUCIF 1 IIIA CZM (I) 420 # 2M S41 )?d 


ZMCZZM( S43) 420 HZM( 42) 
IF €2MC) 69470269 


VRS OZM (43) FUCIF2 0740 420% CUT IF 204 C41 DIZ ZMC 


Ge TO 71 

User, 
ZM(CJ+1]) 82M CI+1 + 2M 0542) 
Ctrl) 3(QCI] 24065942)97200 
CLS] se OCS( S41) 4C5( 423972 
UtJ)sus 

UlCIJ+4]) 3UR 

xO J*L =x +2) 

JIxXSUe2 

DO 72 JCEZIKe FINS 

XC UC) eXCUC+1) 
UeuUCPsULUC4ED 

ZM( ICP EZM( ICO]? 

Vescrev¢ gel) 

Ec scp sE (cel) 

O¢4C)2P( JC4+1) 

QCICVFQC C4) 
CSCICIVSCS(UCH} ) 

ZONT TINUE 

JF ENSUE IN|] 

IF OCNJENIT) 76076073 

DO 7S Tle teNvEnlt 

IF CJMPQJEDITOCEITY) 76e7G4 75 
YENITOCTP VSVEDITCL P= 


ie) 


49 


i Ow 


“€ 


TOT Psa 

JOM AMZ IOVAX =] 
JETARZISTAR]] 
JPL Cs JPL C=} 

Jax Oe IKK | 

DD Fe vps peN’@TOUS 

PE GUT MOM mt) Pee 774 77 
JONTOM EPS JARO pal 
DINT Tue 

IMNDs 2 SLL APS: 

Pe Gey BPM Hd 

PF fy sPRC) BE oe Sle 8} 

Th (J+a—eJOMAK) B24 Rd BG 


ae ode | 

TF (Je SAND OMY) 64972063 
Mat ] 

a0 TO 63 

JOR IE ETNOMUE IN 

SO EN 


PINT LeTEVE eNeyRZeITTAR 
CONT TUE 
IE (SCONCH) 27,870R86 


50 


Parapet ©} 
ee 
-430 
tc JOC 
370 
S289 
8299 
S400 
5410 
3420 
S430 
5440 
Seso 
5460 
3470 
2440 
3499 
S500 
£510 
55e0 
4530 
5540 
Lahet o] 
6860 
5570 


t) 


“ 


¢ 


Fes eres 


EOTT SUBROUTINE FOR THE 2UFF HYDRODYNAMIC CODE 
SUBROUTINE EDIT ‘ 


COMMON C5 (SCC) eCUSP1 (6) eCUSPA(6) eCUSPC 16, sCUSPDI6) eCUSPE(6). 
ICU“PS(6! sOLSCET(I12) 6B (500) eEGSTC(6) sEQGSTDI6) eEQSTE (6) eEQSTGL 
$6) eEGSTH(6) sEGSTNI G6) cEGSTS(G) eJAND(6) eYEDIT(10) eJORG(10) eMATL(I6) oN 
$7 C20) oP 1590) sPMIN(6) 1G SOC} eGHO(6) «RZ620) e TEDIT(25) sULSOG) eVI50006 
OXOET Ie PM CRED: 

COMMON CK OTN OTN UCK UC YES CSF ING UF IN2Z « UPMAXs URZL eUSTAR OUKKO. JZF 
PUL «LIND pi eNJECODT TAN“ TRLS eNPRINSNRZ eNTAPE eNTEDT ePMAXe TIME oe TStUZERDe 


CUT AOD . TPUL 


FORMAT (1LHC6keSHCYCLIBX ee GHT IMEBX «e SHO TNH7X « SHUS TART oe SHUPMAXE Xe GHEMA 
1xK6X oe 6HxX JOMAXSX ¢ THMVOUL SE OX s GHMVEREC aX e BHMOMENTUM/2% 5 110¢2E12 64> 
223100521 204/8X4HDTPREX e6HOTPUL SEX « BHE TOTAL 6X e GHE! NE GAS e BHEMVPCSSX 
24 PHL ~FOUNDS& 6 BHX (CJANDI 24x BHX ( JEND2 5X s THR-BOUNDBK se GHUFINSGE lZede 
4Iley , 

FOOMAT (Li) ) 

rHEINARY DATA STORAGE 

KPA > CIPMAX) 

1F (C#TAPE}) 3¢743 

ISTALDzYCTARe] 

weYETE TAPS GaNe TIME e (OISCOTC To Esl elO) eISTAR se UFING JOM AK Oe ISTARD 

IF CECE «6) Sad 

WEETE TAPE Cel eXl I) UCI POI) oO UI eE CIV eV) CSCI ee leUSTARD) 

IF (SOF 6A) £47 

OPIh:T 6N 

FOAMSAT (164 NEw 06 AT CYCLETIC) 

ENY FILE 6 

DEIN & 

PAUSE LAAGE 

cc TO 3 


MOMENTUM AND ENERGY CALCULATION 


EMVNEG=7 ¢ 
E“vecse te 

FSUMSS 4 

ExSuM=t, 

YJETAQL SUS TAR) 

9 Mt z2eJSTARI 
EVVvsIMes#CUCIIM+UlC VHT IIS! 6 

TF CEMV) Ba Ged 

EMVNE GSE MUNEGCEMY 

ce TO 14 

CMVPOTSEMYPOSS+EMY 

CONT ENGE 

“TTUIMe mW CUM4E (2M (0I)78 0.18687 

EROOGMSE RT SAS PM CIDE CUC II F4UCIO LD EH2 746 LBGE 7/96 
LENTINGE 

eTOTALESSUM+ERK SUM 

Us JO MAW ee] * 
IN VPL=l 

MV EL ET VVEL UC UM) ee Se C 74 SMO] 242MM) ) 
TEC UMe JOVAK? TA0e1Gela 

IF CUCUMel)) 16416014 

IF CUM ep THe b Sede 


hy te 


we 


ome ty 


~ 


Js UMe] 6149 
ce TO 12 61660 
J'!= JPMAX +44 6170 
cMVERE. 6140 
ECuMUYVOO SEVEN 4( UM) (Se RO FMOUM IDO FM (IST DD 5156 
Te CUMeESFIND 1THelE el 6200 
Jus Mel 6710 
to TO 17 6.20 
CUVERZEMVPL eS NVPS oe80 
CTPP2elyPosevax 2 6649 
DIOUL © SEMYPL SEVAY 6250 
42600 

CPINT OUTPUT VARIABLES 6270 

acho 

IONDD2 DOV) 6296 
JONO2Z 2 AND 2) 63500 
PRINT PeNe TIME eO THe US TAs UPMAX ePMAX e KIPMAX AEMVEL cENVERAEMYVP™ TOD 6210 
eNTPULS ee TTOTAL ei MVNEGOEMVPOS e XO1) eX CURNDI 2 eX ( UUND2Z) XC YFIND e JF IN 6320 
LINFZFLINT +5 6230 
TE CXMODEC LIND e FO}? 2lel2oe2l] 6340 
CCImNT > 6350 
RETURN 6260 
END 6376 


$2 


fk eae 


Yr 


ATION OF STATE SUFROUTINE FOR THE PUFF HYDRODYNAMIC COLE 
SUBPOUTINE FOST CE1efiePlem) 

COMMON CS(£00) sCUSPI (6) eCUSPA(6?) sCUSPC (6) sCUSPD(6) eCUSPGI6) 6 
1CUSPS (6) cDIECPT C12) o£ (800) eEGSTCLE) eEOSTDI6) eCQSTE (63 sEOSTCL 
26) eESETHIG) SEGSTN(G) s0EGSTS (6) eJEND(6) eYEDITOCIDO)D eYORGI10) eMATLIO)e & 
BZ625-79F 0500) ePMINGE) 630500) ARHO(G) R7620) eTESIT( ZED eUL SOO) VENI s 
GMOS 03a PE (RTO) 

LOMEON Unbee DT Ne DT NM UCK OIC YC S euF ING UF ING e UPMAX e JRZL ce USTAR INKL 6 2 
LUL sL INTs NeNJEOLTeNMTOLS eNPRINONEZ eNTADE ¢NTEDT ePMAXe TIME  TSsUZCPOs 
2wTsel,7P 


CNG SD CIA 
FMtirOttiele 
VIE™O CMY sD 

TF CES) 1eSeF 


VAPOR EQUATION 


ENUSFOGCTNOMD® (levi) 

TF OCTN Gt 4to 9) Fe Fe? 

TPSTIT tl OM) HO eM XRF CE NUS) ) 

oc TC & 

Tryarr ante (yy) 

TO OSE NUH CLOS TH iM) + CEQSTGIM) -EQSTNHIM) ) *#SOPTF SCENU)) 
PLS4AXTECPMING | Pe CF LTE PHT S2SRRHOCNM) ) 


se Ta 4 


TwoRKwAVE SOLID SQUATION 
IT CCUSMACMD) Bebe 
ARCSTVURSTUSPAIY) 
TFC AGC y We Ge7? 
TOLECOTUSPS OM) #3 RG+CUSPE CS) ) RAORG+CUSPT (IM) ) HARG 
CLPSECUTH PL OMP+T OO SeZt SCUSMSCM) ED 
Oke Fe Oe, 


ae 


DNEWAVE SSLTO CQUATION 


TPE COTA TS OM) HEME OSTLOM) DHEMUFEGSTC OM) ) REMY 
BypelITFel Pel Vato coup ed 
9° TURN . 


Oty 

270 
6409 
aie 
Gal 
Gard 
£400 
Hea 


—~ a 
ID AGE 


i™ 

aati 
yi 
. 


G 
nO 


+ 
2 


sos 
GG 
-@Q 
aD 
* G5 3H ada ETAT, RS PTS SLA LOT TRIS 8 EE RRND BD ee REAM I. NeeT Se 


eee oe 
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APPENDIX ill 


(PUFF } 


Cia ME og One eo eae 2 
82 
; JBNO au 
. fee Tt EE 
5 il 7 
| | 
| 


rr ec Tn a NT 


be See 
4 ANGLE 
3 NRZC NMTRLS JR2U J2PUL.~—s*NPRIN 
3 2 ICO 40 25 
2 OTScPT (I), 271, 10 


MH TESTPROBLEM NGL MRM TWA MATERIAL XH ARBITRARY SPECTRUM . 


TABLE DECK FOR PROBLEM 1 


/ 
NHNU NTEDT NJEDIT L@#HI2 
0) 0 \ 


109 
1 


BA? A} B(28) ENGEi2@) AA(2.5) 8(2,5) EDGE(2,5) AA 


DGE(23) 
4132+] 


?2 6951+ 2 -1.944 4.958+1 T2342 -1.562 619841 279341 
A(2a) Bl ) EDGE(2 1 AAI2 2) B(2 2) FOcE2 2) AA(2, 3) B(2 al 
2l S66 89 1.535+4 re 78 3.09941 3161 3161+ Sees 5, 
SoM CUSPHET CUBE CosE a cupReTay—cusretay CUSPS(2) i 
20 0. 
19 me i) — OSES) a ava) “ QSTH(2) EQSTS(2) PMIN(2}] 
2699 B54 -1232+12 122¢iL 2.04 25 0 6. +8 
Pag Sa 


i” 
| 


: 6(}.9 DGET1,9 Usd B(L10 SETI IG AA(LI : 
16 1.7121 -58565-1 49886441  74648-1 -3.5208-1 6.198+1 5.6627-1 -2.9836- 
Rie BOE A ih, BOEli 7) AALS BTi.5) OSE 
15 -1.5975 24798+1 140214) -2. 1653 3.09941 3.3035 -74228-1 4.1321 


aerels LAC ay ia ee Ss By ee if 
14 / 1239641 225942 -2.0539 1849841 2.25142 ~2.0839 2.06641 4880841 
Bt1,3 


AAGL,L) Bil.}) EDGE( 1 eae (12) B22) E0GHI2 AALS 

13 ( 6565+2 -28569 1870-1 6612+2 -8. aay 2040-i 1.8784 
CUSPI()) CusPa(i)  CusPct) USPDt USP G1 2-16 
2 0. 0. 0. 0. 0. 0. 


RH@(l) EQSTC(L) EQSTD(:) EQSTE(1) EQSTB(1) EQSTH()) EGSTS(1) PMIN( 1) 
22 2.513411 1668411 7559+10 78 25 1.748413 -1.4+9 


of 


UCYCS  NTEST 
25 ¢ 250 30 
CKP a ce 

24 10 2.-6 


93/  -1.121 826443 2581 -SA76-l1 150+2 


— 


58 


B26 Ge 6) RATA TY B27) EDGEIO TY 


TALL DUCK FOR TLET FRODLEM i 


2L1567. TORMS 


TBLEI) LI(1) THL CARL # 
0.5972897L 00 0, l 
0.7712318L 00 0. 2 
O. 8991189 O60 ihe 3 
0,1004823E 01 0, 4 
0.10970181 o1 O, 5 
G.1179999E Ol 0. 6 
0,1256234E 01 0, 7 
0.1327283E 01 0. 8 
0,1394228E 01 0. 9 
0.14S7829E Ol 0. 10 
G,1518656E 01 0. 11 
0,1577137E 01 0. 12 
0.1633631£ 01 0. 13 
0.1S8B417E Ol 0. 14 
O.1741701E 01 0. 1s 
0.1793693r 01 0. 16 
O.1844536E OL o, 17 
.1894372L 01 U.5026790E-30 18 
U.194330SE 01 0,.945556SE-28 19 
0.1991455C 01 Q,10184€7E-25 20 
0.2038300E 01 0.6780110E~24 21 
0.2085728L 02 0.3002525E-22 22 
0.2132004E 01 0.93258S6E-21 23 
0.2177796L 01 0.2127423E~19 2i 
0.2223151 ol 0, 3701496E-18 28 
0.2268135i 02 0.5078059F-17 2 
0.2312798L 03 0. 564S669E~16 27 
0.2357167E 01 0.5202916E-15 28 
0.2401296L 02 0.4058228i-14 29 
0. 2445216L Ul 0.2725063E-13 30 
0,2488976E 01 0.1600093E-12 31 
0.253259¥E 01 0.8319149E-12 32 
0,2576116E Ol 0.3876429E~11 33 
0.2619573E 01 0.1635183E~10 34 
0.2662981E 01 C.6297824E=10 35 
0.2706380E 01 0.2233014E=09 36 
0.2749789E 01 0.7338907E-09 37 
0.2793255E Ol 0.2250884E-08 38 
0.2836786L 01 0.6475 984E~08 39 
0.28804130 01 0.1756749E-07 40 


Cee eta 
Os 2308066. 
0, 307271 28:, 
U, 3UDG4 tun 
O,310092L 


yA ef OL 
Ure va ea a, 


O.31907002: 
0. 32 36546L 
Q,32817140 
0.3327763L 


0.3374200L 


O.3421056L 
9. 34E8363L 
0. 3516174. 


O.35645038L t 
U.3613382L i 


0. 3628537 
Uadadeo rel 
U.376377 21. 


OC soloed TEL. 


D.3KE754h, 


Und 206I3E. 


De STH 5aN. 
OeUZINSHL 
O.4065 3905. 
QB P42 31a: 
O.4200 S84l 
42596771: 
O.8320261L 
0.4382229F 
O.4445692L 


O.GS1L0757L 8 


O.4S77551E 
UH6NE 9 3L 


O.4736850L 


O.4789699. 
C.4 8649161 
U.G942 7141. 


G. 5023761. 
O.SLO7L14b | 
0,51943131, 


9,5285303L 
O,5380537L 


CG oesuouaL | 


9.558883 73, 


Co 
r4 
o 
4 
we 


. 
a“ 
c 
rer 


mot 
e 
ar fe or 


Dh ae Gd tO 


iG 
co cd 
me 
woe is 


(2 
. 


Ow Orws Mee 
C 
an 


47 8D ps Ut pe 
wan Nd SS 


a) 
. 


0.33704330-04 
0.599232 aL-048 
0.1637839L-03 
O.1755874r-0¢4 
0.290€129::-03 
U4711179L-03 
OQ, 74890891'-03 
O,1168785:-02 
O.1 792640 -02 
QO. 27048741,-02 
O.4017908:.-U2 
D.S8$8L003-02 
U,R4R9ESTE}07 
O.12604 54-01 
O.170171451.-01 
QO. 2266081 -01 
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28 May 1965 
ADDEN DUM 
AFWL TR-65-248 THE PUFF AND P PUFF COMPUTER PROGRAMS 


This addendum is to notify users of recent improvements and changes to 


the PUFF and P PUFF programs and to correct two minor discrepancies in 
AFWL TR-65-24. 


The improvements are related mostly to REZONE and should help to improve 
pulse resolution in poorly zoned problems. The major changes are found in 
the Equation of State Subroutine. The solid phase equation of state is now 
fit to raw Hugoniot data. 


These additions will affect the solutions to the three sample problems 
Slightly. If one wishes to check the code conversion to his computer, we 


suggest adding the following changes after the initial conversion has been 
checked. 


On page 23, delete the firs* statement past statement number 17, SSCK=1.. 


On page 2h,delete statement number 53 and substitute the following 
three statements: 


GO TO 531 


53 SSCK=1, 
531 DIN=DINH . 


On page 31, change the first statement before statment number 2, 
IF( JPMAX-J ZPUL~J RZL-JPMAXI) 8&7, 87, 2 
to read 
IF( JPMAX-JZPUL-JPMAXI) 87, 87, 2. 


On page 33 delete statement number 50 and one statement past, JRZV=1 
and MRZ#1. Number the next statement, JFINO=JFIN, 50; 50 JFINO=J FIN. 


Change the second statement on page 34 to Mel. 


On page 34, delete the first statement past statement number 52 and 
substitute the following four statements: 


IF(TIME-SDUR) 520, 530, 530 
520 KRZ=) 

GO TO 53 
“30 IF(JPLC-JV-JRZL) 53, 53, 60. 


On page 34, change the first statement before statement number 56 to 
read 


LF(JPLC+1-JBND( MM) ) 57, 56, 56 - 


On page 35, delete the second statement past statement number 86, NRZaN, 
and substitute the following three statements: 


IF( JPMAXI) 860, 860, 870 
860 JPMAXI=0 
870 NRZ=N . 
On page 38, change statement number § and one atatement past to read 


7 TS2=(CUSP1(M)+( ( CUSPS(M) ®ARG+ CUSPD(M) ) SARG 
+CUSPC( M) ) “ARG) #(1.-(CUSPG(M)® EMU) /2.) 
Pl=TS2+CUSPG(M) *D®E1, 
and change statement number 8 to read 


8 TS2 = ((EQSTS(M) #EMU+EQSTD( M) ) FEMUFEQSTC( M) ) 
#EMU"(1,-(EQSTG(M)*EMU)/2.) - 


On page 41, delete statement number 20 and substitute the following 
three statements: 


20 IF(N-1) 200, 200, 210 
200 JPMAXI=#JFMAX+10 
210 IF(J-JSTAR) 22, 21, 21 - 


On page 44, change the second statement past stetement number 31, 
NRZ=-50, to read NRZ*50. 


On page 46, change the first statement before statement 1. .mer 2 to 
read 


IF( J@MAX-JPMAXI) 85, 85, 2 - 
On page 48, change the first statment before statement 54 to read 
IF(JPLC+i-JBND(MM))55, 54, Su - 


On page 50, delete the first statement past statement number 84, NRZ=N, 
and substitute the following four statements: 


JPMAXI®JPMAXI -J RZ 
IF(JPMAXI) 8b0, S40, &50 
BLO JPMAXI#0 
B50 NRZ=N . 


On page 53, delete the first executable statement, ENU=D/RHO(M), and 
substitute the following five statements: 


IF(EL) 10, 10, 20 
10 E10, 

P1sd, 

GO TO 9 
20 ENU=D/RHO(M) . 


Also, there are two more changes on page 53. Statement number 7 and one 
statment past are changed to 


7 TS2=({CUSP1(M)+(( CUSPS(M) *ARG+CUSPD( M) ) *ARG 
+CUSPC({ M) } ARS) #(1.-(cUSPG(M)® EMU) /2.) 
P1=TS2+CuSPAa( M) ®D*El , 


and statement number 8 is changed to 


8 TS2=( (EQSTS(M) *EMU+EQSTD(M) ) “EMU+ECSTC! M) ) MEMUM(1.-(EQSTG(M)*EMU)/2.). 


This addendum has been reviewed and is approvec. 
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